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£h '• Abstract 

ci 

We show that an explicit counting of Gribov copies can shed light on the infrared 
behavior of non-abelian gauge theories. A power-law growth of the number of copies 
suppresses gluon propagation while the distribution of copies along a gauge orbit implies 
an enhanced density of small eigenvalues of the Fadeev-Popov operator. Both of these 
phenomena are related to confinement. The discreteness in the number of copies and 



in 



Or 



■ the associated nonlocality also has implications for vacuum energy. 

i> 1 Introduction 

o\ 

q ■ The continuum definition of non-abelian gauge theories is sensitive to the global structure 
of the field configuration space, since gauge equivalent configurations remain after standard 
^ ; gauge fixing procedures are applied [I]. Rather than just a technical nuisance, the structure 
of configuration space that is implied by the existence of Gribov copies may be central to the 
understanding of a mass gap and the mechanism of confinement in these theories. 

A gauge fixing condition G(A) = defines a hypersurface in the gauge field configuration 
space. But a gauge orbit, a set of gauge equivalent configurations, can intersect this hypersur- 
face more than once. To characterize this we consider the following integration along a gauge 
orbit where A v n = WAM + ^WdM, 
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1 + N(A) = j [dU}6(G(A u )) det (^P) 



(1) 



This result has an analog in ordinary calculus. 5G(A U ) / SU is the Fadeev-Popov (FP) operator 
which we henceforth denote by T>a- N(A) is the number of Gribov copies, the number of 
additional intersections of the orbit with the hypersurface, and it is clear that N(A U ) = N(A). 
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One source of Gribov copies arises from the existence of different regions of configuration 
space in which the spectrum of T>a has different numbers of negative modes. At a boundary 
between two such regions an eigenvalue of T> A is changing sign. Gribov proved that pairs of 
copies can exist close to boundaries, with one copy on each side. The first Gribov region T>a 
has no negative eigenvalues, and this region includes the perturbative configurations. Gribov 
copies can be ignored in the perturbative limit. 

The generating functional for the quantum theory can account for the presence of Gribov 
copies as follows p], 

Z[J] = j [dA]e iS[A ' J] 5(G(A)) \det(V A )\ { - 1 . (2) 

S[A, J] is the standard pure gauge theory action with source J. Given that "this problem 
seems to be almost hopeless" jl] given the presence of N(A), attention was focused on attempts 
to directly excise copies from the functional integral. Of particular interest is the fundamental 
modular region (FMR), which is a subregion of the G(A) = hypersurface such that every 
gauge orbit intersects it once and only once. It is contained within the first Gribov region 
although there is some overlap of the boundaries of the two regions [3]. 

From the definition of the FMR we can write [dA] = [dA FMR ] [dU] (see jl] for more details), 
and then from (CQ) we have 

Z[J] = j [dA FMR ]e lS ^ J \ (3) 

This restriction on the functional integration is difficult to implement in a continuum de- 
scription due to the highly nontrivial nature of the FMR. Instead a restriction to the larger 
first Gribov region is often considered since it can be implemented as a positive definite con- 
straint on det (T>a)- Then the absolute value in (T5]) may be removed and it is also argued that 
remaining copies are such that the 1/(1 + N) factor can be neglected [5]. 

Our approach in this work will be to remain true to (j2J) and allow the functional integral to 
sample the Gribov copies. A prototype of this type of approach to Gribov copies was studied 
in detail in [6]. One goal then is to obtain information about N(A) directly through explicitly 
counting copies. In addition we shall argue that the distribution of copies along a gauge orbit 
provides information about the spectrum of T>a- The behaviors of N(A) and Da that we find 
have implications for the confinement mechanism. 

The Landau gauge condition G(A) = d^A^ = gives V A = —d^D^ with the standard 
gauge covariant derivative while the Coulomb gauge condition G(A) = V ■ A = gives 
T*a = —V • D. Gribov copies exist in both gauges. Coulomb gauge lacks Lorentz invariance 
and is in fact is not a complete gauge fixing even without Gribov copies. But it provides a 
more direct connection to the confinement mechanism through the long distance behavior of 
the instantaneous color Coulomb potential, where the latter is governed by the low momentum 
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spectrum of T>a- Also, since ghosts are absent in the Coulomb gaugej 1 ] we can infer that the 
ghosts present in the Landau gauge will also be absent from the physical states. This must 
be the case even though (T5]) lacks a Becchi-Rouet-Stora-Tyutin (BRST) symmetry. In fact it 
is presently not known how to maintain a local BRST invariance while accounting for Gribov 
copies jH |7|. BRST symmetry emerges in the perturbative limit where Gribov copies are 
absent, and then the construction of the physical states can be accomplished explicitly. But 
the decoupling of ghosts is more general. 

In addition to being gauge invariant, N(A) is also scale invariant. For a solution to 
d^A 1 ^ = for example there is another solution on another gauge orbit where A(x) — > XA(Xx) 
and U(x) — > U(Xx). Then the two gauge orbits are related by a scale transformation and the 
Gribov copies on the two orbits are in one-to-one correspondence. Thus gauge inequivalent 
configurations related by a scale transformation have the same N(A). 

An explicit counting of copies is made possible by a restriction to configurations of finite 
norm in Euclidean space, 



Every gauge orbit does have configurations of finite norm and so this constitutes a partial gauge 
fixing. The allowed gauge orbits are now characterized by the value of the norm functional 
and at the stationary points the remaining gauge fixing d^A^ = is satisfied. This makes 
the number of copies N(A) well defined. The norm along a gauge orbit has reached a global 
minimum at the point it intersects the fundamental modular region by definition. In other 
words the gauge fixed hypersurface restricted to finite norm contains within it the FMR. Our 
focus then is on the 1/(1 + N(A)) factor that is needed to compensate for the copies in this 
space of finite norm configurations. 

In particular we would like to determine how N(A) depends on the amplitudes A^(k) and 
wave vectors k of the fourier modes in a finite Euclidean volume V = L 4 , 



Note that once the volume is fixed then we can no longer consider scale transformations, since 
the spatial extent of a finite norm configuration should also change under a scale transfor- 
mation. In section 3 we study configurations that are approximations to finite volume plane 
waves and which allow a counting of Gribov copies to be performed. In section 4 we study 
the distribution of copies along the gauge orbits and discuss the connection to confinement. 

1 The cancellation of | det(2?^)| due to the Aq integration is described for example in [S]. 
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2 Implications 



In this section we shall discuss the implications of two of the results that follow from the 
counting of Gribov copies, as described in the next section. 

• For A^tyA^l— k)/k 2 less than some constant, N identically vanishes. 

• When N is large it grows approximately as a positive power of A M (fc)y4 A '(— k)/k 2 . 

The first result is a consequence of the discrete nature of N and it means that there is a 
region within the space of finite norm configurations for which N vanishes. The gauge fixed 
hypersurface within this subregion will be a subregion of the fundamental modular region. 

Given these properties, we would like to illustrate the effect that the 1/(1 + N) factor in 
the functional integral can have on the gauge-fixed propagator. For this we shall develop a 
toy model that shows the effect on the propagator of an otherwise free theory. For a free 
theory the generating functional in momentum space can be written as a product of simple 
one-dimensional integrals, one for each k. We replace A^ik) by A and consider the Euclidean 
form of the integrals where a space-like momentum is mapped to a positive k 2 . And further 
we shall temporarily suppose that the 1/(1 + N) factor also factorizes over the momentum 
modes; we shall comment more on this below. We thus consider the toy generating functional 
based a one- dimensional integral: 

f°° 1 1 
Z(j)= dA u exp( — k 2 A 2 +jA). (6) 



The toy propagator is 



j=0 K 



The dependence of iV on ^A 2 follows from the properties of N noted above and dimen- 
sional analysis. We shall model this dependence as follows: 

A 2 <£ 

The appearance of a physical mass scale A corresponds to an explicit breaking of scale in- 
variance; it is a physical scale below which Gribov copies are important and above which 
they are not. From (Q we see that A determines the cross over between the infrared and 
ultraviolet regimes as follows: for k 2 ^> A 2 the exponential factor constrains the range of A 
that contributes to the integral, while for k 2 A 2 it is the 1/(1 + N) factor that constrains 
the range of A. 
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It is interesting that a mass scale A enters as soon as the classical structure of the gauge 
configuration space is put into a quantum mechanical context, due to the A 2 /k 2 dependence 
of the former versus the k 2 A 2 dependence of the latter. It is the classical effects of the 
gauge field interactions that are responsible for a nontrivial N. The direct quantum effects of 
interactions are being ignored here, but they also give rise to a mass scale through dimensional 
transformation. It seems reasonable to assume that these two mass scales are similar. 

For this model G(k 2 ) can be determined analytically, but we don't give the messy expres- 
sion here. The leading behavior at k 2 <C A 2 is 



G{k 2 ) 



i i 

Fln(A 2 /fc 2 ) 

§ln(A 2 V 2 )^ 

1 q-l fc 2 
3 a-3 A 4 



a 
a 
a 



1 

2 
3 



> + 



a > 3 



(9) 



Our explicit counting results of the next section show that a is certainly greater than 3 and 
is closer to 5. The leading behavior at k 2 3> A 2 is 
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2 A 2 

7T k 2 




(10) 



We can also consider a smooth analytic continuation to negative k 2 , and then we find that 
G{k 2 ) -> (1 - a)/k 2 for -k 2 > A 2 . 

We can compare these results to the original analysis of Gribov [TJ [9] , which described an 
attempt to remove the 1/(1 + N) factor by instead restricting the configuration space to the 
first Gribov region. The implementation of this restriction in his semi-perturbative approach 
had the effect of introducing into the functional integral another factor instead of 1/(1 + N). 
This factor happens to factorize over the fourier modes and in the notation of our toy model 
takes the form exp(— pr4r). This produces the propagator 



G{k 2 



A 4 



+ 



k 2 + W 



4 (i - K + 
p I fc 4 + ••• 



P<A 2 
k 2 > A 2 



This propagator vanishes with k 2 in the infrared, which is also the case for our model for 
a > 3. Although this infrared suppression of propagation is different from a conventional 
mass, there are no longer massless excitations in the infrared and in this sense a mass gap has 
formed. The vanishing of the propagator implies that it violates of positivity [10], and in the 
context of gluon propagation in QCD this is expected for confinement. Violation of positivity 
is also found in lattice studies [TT] . 

The main difference in the two models is in the ultraviolet behavior, where our propagator 
approaches the free propagator exponentially quickly according to (flQl) . in contrast to the 
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power law approach in (JTTJ). A power law approach to perturbative behavior is common to 
all versions of the Gribov-Zwanziger framework [TJ, [121 UHl [7] . 

An exponential approach to the free propagator has a significant implication. The operator 
product expansion for QCD with massless quarks relates the asymptotic behavior of the gluon 
propagator to various condensates. 

/, 2 net 2^-1 fc2 ^°° , {G^vG^) (qq) 2 

(k G(k )) ~ ai + a 2 — ^ ha 3 ^- + - (12) 

Thus if a gluon propagator has ultraviolet behavior as in ( flOl) or (JTTJ) it is associated with a 
vanishing or non- vanishing gluon condensate respectively. The gluon condensate in QCD is 
in turn related by the trace anomaly to j(Tjf), which by Lorentz invariance is the vacuum 
energy. Our model yields vanishing vacuum energy, and at a superficial level this could be 
said to reflect the classically scale invariant nature of Gribov copies. 

We have traced how the quantity A that reflects the breaking of scale invariance enters the 
description of Gribov copies. One might have expected the breakdown of scale invariance to 
be manifested by a non- vanishing gluon condensate of order A 4 . But we see that this is not 
the case. In fact we can relate the absence of a gluon condensate to the discrete nature of N; 
because A" is discrete it can and does remain identically zero until ^A 2 is sufficiently large as 
in (JHJ). This is the cause of the exponentially small effects of Gribov copies as in (fit)]) in the 
regime k 2 ^> A 2 since the contributions from A 2 ^> 1/k 2 are exponentially suppressed. The 
effects of Gribov copies are thus not sufficient to produce a local gluon condensate, which as 
the operator product expansion demonstrates is a modification of the short distance behavior 
of the theory. This reflects the nonlocal and long distance nature of Gribov copies. 

Finally we comment further on the likely false assumption that 1/(1 + N) factorizes in the 
space of fourier modes, which led to our toy model. For the toy model the region over which 
A" identically vanishes is a hypercube in the space of all A^k). In a more realistic description 
the hypercube is probably replaced by some other shape. But it clear that the reasoning of the 
previous paragraph depends on the mere existence of a region where N identically vanishes, 
and not on its shape. Note also that in QCD the coupling falls in the ultraviolet, and thus 
the effect of Gribov copies in the ultraviolet is even less than what our toy model indicates. 

3 Counting 

For our purposes here we shall consider static gauge field configurations, in which case the 
Landau and Coulomb gauge conditions become identical. We restrict further to static spher- 
ically symmetric configurations. As we shall see, these configurations have finite A" when 
restricted to finite spatial norm. We take the gauge group to be SU(2) with gauge field 



6 



A^(x) = A®(x)a a /2i. Gauge configurations of the form 

A?(x) = (13) 

automatically satisfy the gauge condition diA\ = 0. By an appropriate choice of f(r) we 
shall be able to approximate plane waves having only spatial oscillations, and thus model the 
spacelike k 2 of the previous section. 

The following gauge transformations produce more general spherically symmetric gauge 
configurations, 

U(x) = cos(a(r)/2)l + isin(a(r)/2)—a a . (14) 

r 

For the transformed gauge field to be a Gribov copy it must satisfy di(A u )1 = 0, and this 
leads to the Gribov pendulum equation [U [91 fT3] : 

rV'(r) + 2ra'(r) - %/(?") + 1) sin(a(r)) = 0. (15) 

The classical scale invariance corresponds to solutions related by f(r) — > f(Xr) and a(r) — > 
ot(Ar). Here we only study the classical configuration space and thus we treat the coupling g 
as a constant; also the mass scale A will not appear in this section or the next. 

From ( [TBI) . f(r) — > as r — > at least as fast as r to have a regular configuration. To 
maintain regularity under the gauge transformation, a(r) — > 2nn for integer n as r — > so 
that U — > ±1. For a given a(0) the set of solutions of the pendulum equation (the set of 
copies) is parameterized by the "initial condition" o/(0). Since the equation is invariant under 
both the constant shift a(r) = a(r) + 2mr and the reflection a(r) — > — a(r), and since the 
g dependence can be absorbed into f(r), without lack of generality we may set a(0) = 0, 
a'(0) > and g = I. 

The infinite number of copies parameterized by a'(0) reduce down to a finite number, 
possibly zero, under the further restriction to finite norm configurations. For the configuration 
in ( [TBI the requirement is 

//*oo 
d 3 xA$(x)A?(x) = 8tt J drf(r) 2 < oo. (16) 

For this to be also true of H^H, the gauge transformation must have the property that 
a(r) — > 2nir as r — >• oo for some integer n. Since ct(0) = this implies that the gauge 
transformation has a winding number n. The pendulum equation is such that these finite 
norm solutions only exist for a finite set of a'(0). On the other hand the infinite set of infinite 
norm solutions are characterized by a{r) — > (2n + 1)tc as r — > oo [9]. 

Since we work in the space of finite norm configurations, by Gribov copies we shall hence- 
forth refer only to the finite norm solutions. Because (ITS"]) can be efficiently integrated as 
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an initial value problem, copies can be counted by simply observing the behavior of a(oo) 
as a function of a'(0)§| For example for two values o/(0), a(oo) may jump from some value 
(2ni + 1)tt to another value (2ri2 + 1)tt. We then know by continuity that there are at least 
| ni — rz.2 1 copies between these two values of a'(0). The numerical task then is to calculate a 
set of a(oo) for a sufficiently fine grained set of a'(0), and then obtain a value for N by adding 
the | ni — rial for all the jumps in a(oo). Copies can be missed when a(oo) jumps forward 
and backward within the a'(0) stepsize. As N becomes larger the copies are becoming closer 
together, and thus the stepsize of a'(0) must be appropriately small to approach the true N. 
In addition the maximum value of o/(0) beyond which there are no further copies is increasing 
for larger N. 

We shall use these spherical configurations to approximate spatially dependent plane waves. 
We thus consider configurations that oscillate as a function of r and which are modulated by 
an amplitude which vanishes in the small and large r limits. When the radius of the amplitude 
peak is much larger than the wavelength of the oscillation we are approaching plane waves. 
For the radial profile we choose the form 



f(r) = Ae r/R (l - e r/R ) p sm(kr), p = 1 or 2, 



(17) 



where we need to choose R large compared to 1/k. k and A are analogous to the k and A 
of the toy model in the previous section. Holding R and p fixed is analogous to holding the 
volume fixed, while a scale transformation f(r) — > /(Ar) is equivalent to k — > Xk, A — > A, 
and R -> R/X. 
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Figure 1: a(oo)/27r versus arcsinh(o; / (0)) for A = 100, k = 1, R = 30, p = 1. 
2 Numerically we integrate a(r) out to r = 10 10 . 
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We choose the discretization a'(0) n = sinh(nA) for integers n and small A, since then 
the copies are spread out quite evenly as a function of nA (or arcsinh(a'(0))). In Fig. (1) we 
present an example of the behavior of a(oo)/27r for A = 100, k = 1, R = 30 and p — 1, where 
we obtain N = 80. Copies (of finite norm) exist where a(oo)/27r crosses integer values. These 
values are the winding numbers of the gauge transformations, and we see in this example that 
there are 10 copies of zero winding number. 

Our counting results show that N is zero unless A/(k + c\) > C2, where the constants 
depend on our modeling of plane waves. We plot values of A versus k at which the first copies 
occur in Fig. (2) for (R, p) = (30, 1), (100, 1), (30,2). The linear relationship holds very well 
and we find that C\ ~ 0.45,0.18,0.0 respectively. We understand the offset C\ as follows. For 
small k we are not approaching a constant gauge field, but rather a spherical configuration 
with a characteristic size R, corresponding to a superposition of plane wave modes. c\ tends 
to zero when the radial profile is more strongly suppressed near r = 0, as when p = 2. 




1 L5 2 2^5 3 3J 4 



Figure 2: A versus k at which the first copies appear for different models of plane waves. 
{R, p) = (30, 2), (100, 1), (30, 1) from top to bottom. 

As motivated by the last section, the main point of this section is to find the behav- 
ior of N when it is large. To investigate this we study the dependence of iV on A for 
large N and fixed k. In Fig. (3) we display log(iV) versus log(A) for the cases (k,R,p) = 
(2, 30, 1), (1, 100, 1), (1, 30, 2). The number of copies we have calculated per configuration 
ranges up to nearly 14000. To achieve a single value of N this large we generate over 350000 
numerical solutions of the Gribov pendulum equation, § The roughly linear behavior seen in 



3 We used the rk8pd method of the GSL C package. 
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Fig. (3) indicates a power law growth iV oc A a with a m 4.9,4.9,5.5 respectively. This power 
law behavior leads to the model we investigated in the previous section. 




5.2 5.4 5^6 5^8 6 6.2 6.4 6.6 6^8 7 ~1~.2 



Figure 3: log(iV) versus log(A). (k,R,p) = (2,30,1), (1,100,1), (1, 30,2) from left to right. 



4 Towards confinement 

Since the Gribov pendulum equation has a symmetry a{r) — »■ —a(r\ every copy is duplicated 
under a'(0) — > — a'(0), and thus the total number of copies is eveno Let us continue to focus 
on the copies with a'(0) > 0. Consider the values of a(oo) for a'(0) — > 0+ and a'(0) — > +oo 
respectively and let us refer to these two values as «o(oo) and a^oo). For example in Fig. (1) 
both values are +7r. From our numerical analysis we find that 000(00) is always +ir while 
«o(oo) can be ±7r. When there are no copies ao(oo) = ir. Often there is a range of A for which 
only one copy exists (along with its duplicate), and this copy first appears (as A increases) 
when «o(oo) changes from 7r to — tt. The copy has a value of c/(0) which is infinitesimally 
close to zero for some value of A. 

In this case the field configuration is related to either of its close copies by an infinitesimal 
gauge transformation U(x) = 1 + X(x) and this implies a vanishing eigenvalue of the FP 
operator, T>a = —diDi, 

8 t Ai = = dAU^AiU + -U~ l dm = --V A X(x). (18) 

9 9 

4 The TV's given in the previous section ignored this duplication. 
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In such a case the curvature of the norm functional is also vanishing, 



\A\\ + 2 J (PxtiiXdiAi) + J d 3 xti(X j V A X) + C(X 3 
\A\\ + C(X 3 ). 



(19) 



Thus as A increases the first Gribov boundary is being crossed and an originally positive 
eigenvalue T>a is turning negative. An analogy is to the function (x 2 + a 2 ) 2 when a 2 changes 
from positive to negative; the curvature at x = turns negative and two new minima ap- 
pear. Thus we expect the original configuration to pick up a negative eigenvalue while the 
two copies do not, and so the copies remain in the first Gribov region while the original con- 
figuration moves outside. The two copies also remain in the FMR, which thus incorporates 
a reflection symmetry, and it can be checked numerically that the copies have smaller norm 
as expected. The crossing is occurring where the FMR and first Gribov boundaries coincide. 
This bifurcation picture is described in [2]. 

As A increases further Gribov copies can continue to appear infinitesimally close to the 
original configuration, and thus more Gribov boundaries are being crossed and T> A gains more 
negative eigenvalues. Thus increasing A is pushing the configuration into ever higher Gribov 
regions. For the case R = 30, p = 2, k = 1 we plot «o(oo) versus A in Fig. (4), which thus 
shows the values of A at which Gribov horizons are being crossed. 

,5j 



1.6 



2.2 



-5H 



2.4 



2.6 



3.2 



Figure 4: a(oo)/27r versus log 10 (A) for qj'(0) — ► 0+. 

The gauge orbit can cross a Gribov boundary anywhere all along the orbit, and this 
corresponds to a pair of copies appearing at a finite a'(0) (along with a duplicate pair at 
— a'(0)). The difference in a'(0) for the pair starts from zero and increases as A is increased. 
The analog now is the function x(x 2 + a 2 ) as a 2 turns negative, so that a new minimum and 
maximum appear. This is how most copies appear and in fact it sometimes happens that the 
first copies appear in this way rather than as described above [2]. In this case the original 
configuration stays within the FMR and the new pair is near a Gribov boundary that does 
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not coincide with the FMR boundary. The new pairs may also be related to the original 
configuration by gauge transformations of nonzero winding number. 

By the time there are many copies then it will often happen that two copies are separated 
by a small but finite gauge transformation. Then the X(x) appearing in (fl9l) is determined 
by two solutions cti(r) and a 2 (r) with 0^(0) ~ «' 2 (0) and aii(oo) = a 2 (oo) = 27m. The 
gauge transformation determined by a(r) = (1 — s)a 1 (r) + s«2( r ) with s changing from 
to 1 moves along a gauge orbit connecting the two Gribov copies. This path in configura- 
tion space moves off the gauge fixed hypersurface but it remains within the space of finite 
norm. Thus while X(x) will no longer represent a zero mode of T>a, if it is small it can still 
represent a direction at the stationary point of the norm functional for which the curvature 
C = J d 3 x ti(X^T>AX)/ J d 3 x tr(XtX) is small, \C\ <C 1. The two nearby copies typically 
have curvatures of opposite sign corresponding to being on either side of a Gribov horizon; an 
eigenvalue of Da is changing sign as the Gribov horizon is crossed. The difference Aa'(O) is a 
measure of how far apart the two copies are. We expect a linear relationship \C\ oc Aa'(O) for 
small separations. To illustrate, the representative function x(x 2 — 3a 2 )/6 has two stationary 
points at x = ±a with respective curvatures also equal to ±a. 




Figure 5: a(oo)/27r versus arcsinh(a'(0)) for a large N case. 

Thus if we obtain the distribution of Aa'(O) values for the "close" pairs of copies, on a 
gauge orbit with many copies, we are obtaining a distribution of the small curvatures, which 
are the small diagonal matrix elements of V 'a- We expect these matrix elements to approach 
the eigenvalues in the small curvature limit, and thus we are also obtaining a distribution of 
the magnitudes of the small eigenvalues, the |Aj|'s. This distribution enters the functional 
integral via the | det(T>A))\ factor in ([2]), which is also being sampled at all the gauge copies 
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along the orbit. 

We need only consider pairs of copies having the same winding number, since only then 
can the relative gauge transformation become small. We consider two ways to construct the 
distribution. In the first method we consider the separations of all the nearest neighbor pairs 
of copies. In the second method, for each copy we take the separation to the neighboring copy 
which is the closest. This may count the same separation twice, but this corresponds to trying 
to find the smallest curvature at each copy. 

We give in Fig. (5) an example of a(oo)/27r versus arcsinh(a'(0)) (where the latter is 
evaluated for discrete values nA) for A = 1300, k = 1, R = 30, p = 2, A = 1/14000. The 
resulting estimate for N is 13856 copies (ignoring the duplicates). Each separation is some 
multiple of A. To study the distribution in these separations we plot the log of the number 
that each separation occurs versus the log of the separation, for each of the two methods. 
The results in Fig. (6) and (7) have lines added for comparison with slopes of —1 and —4/3 
respectively. The first method suggests that the density of eigenvalues scales as |A| _1 . We will 
consider this case in the following, although the second method could indicate that there is 
more infrared enhancement. We compare this density to the density of eigenvalues of the free 
FP operator (—V 2 ) which scales as A _1//2 |f This suggests that the density of small eigenvalues 
of T>a is enhanced in the infrared, in relation to the free case, by a factor »s |A| -1 / 2 . 




Figure 6: (log of the frequency of each separation) versus (log of the separation) by method 
1. The line has slope —1. 



5 In the free case dk = d\/(2yX) and the eigenvalues are evenly spaced in k for planes waves of fixed 
direction, which is what we are modeling. 
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Figure 7: (log of the frequency of each separation) versus (log of the separation) by method 
2. The line of slope —4/3. 



For a non-abelian gauge theory, the Coulomb gauge Fadeev-Popov operator T>a — — V • D 
determines the classical instantaneous Coulomb potential |14j . 



The color charge can have quark and gauge field contributions, a h = a^ u — f bcd A c ■ E d , with 
V • A a = V • E a = 0. 

Then the Coulomb self-energy of a static point color charge in a finite volume is propor- 
tional to [TB] 



where F n = (A n |(— V 2 )|A„) is a diagonal matrix element in the Faddeev-Popov eigenstates. 
In the large volume limit with an ensemble average over the gauge field we can write the 
self-energy as [15] 



We view this as an integral over |A| and where p(|A|) is the appropriately normalized density 
of eigenvalues. There is the usual ultraviolet divergence cutoff by X max , but of interest here 
is a possible infrared divergence which is related to the growth of the Coulomb potential 
at large distance. This infrared divergence has been argued to be a necessary condition for 
confinement [15] . In the free field limit where the eigenstates are plane waves in 3 dimensions, 
then A = k • k, p(A) oc A 1 / 2 (from d 3 k oc dXy/X), F(X) = A, and there is no such divergence. 




(20) 




(21) 




(22) 



14 



Our previous results indicate that the density of eigenvalues are skewed towards the infrared 
by a power law enhancement, implying that the eigenvalue density in (1221) is enhanced at small 
A, p(A) — > p(A)A -1 / 2 — > constant. We can also expect an infrared enhancement of F(\). 
The eigenstates associated with the small eigenvalues are functions that are close to the small 
gauge transformations X(x), and we find that these latter functions are typically quite rapidly 
varying with significant overlap with the f(r) profile. Thus the eigenstates typically contain 
Fourier modes with k • k ^> A, and this will enhance (A|(— V 2 )|A) relative to A. 

We obtain a simple model of F(\) as follows. Due to the enhanced degeneracy of the 
eigenvalues of T>a relative to (—V 2 ) at small A, the eigenvectors with eigenvalues of (— V 2 ) 
in some range d(k 2 ) correspond to eigenvalues of D & in a smaller range dX, where d(k 2 ) oc 
(dA)A- 1 / 2 . Then we have (A;|(-V 2 )|A;) = k 2 = jf d{k' 2 ) oc ^(d\')\'- l l 2 , and thus F{\) oc 
A 1 / 2 . This an infrared enhancement of A -1 / 2 relative to the free field behavior -F(A) = A. 

The enhancement of both p(A) and F(\) by A -1 / 2 produces a total infrared enhancement 
factor of 1/ A in the integrand of the self energy ffl2l) , and this is more than sufficient to produce 
an infrared divergence. An enhancement factor of 1/A implies changing from a Coulomb 
potential (1/A: 2 behavior) to a confining linear potential (1/fc 4 behavior). 

5 Beyond Coulomb gauge 

In Coulomb gauge the gauge fixing on one 3-dimensional time slice is independent of the 
gauge fixing on another time slice and so results in this gauge apply equally well when the 
time direction has finite extent, i.e. at finite temperature. The infinite Coulomb self energy 
thus persists at high temperature, even above the expected de-confining phase transition, and 
this is one way to see that the infinite Coulomb self-energy is a necessary but not sufficient 
condition for confinement [16]. This then leads us to consider what happens for a covariant 
gauge fixing such as Landau gauge. 

Our counting of Gribov copies of static field configurations applies equally to both Landau 
and Coulomb gauge conditions. We also recall that V>~^ in Landau gauge is the ghost propa- 
gator in the gauge field background. Thus the quantity [D~^ {—V 2 )V^] XX that we considered 
in Coulomb gauge is the 3-dimensional analog of the 4-dimensional ghost loop contribution 
to a 2-point function, where each vertex corresponds to the ghost current j^ h = c a d ll c a . The 
2-point function in momentum space is 

n^(g 2 ) = - J ^d- a \v)V a \ v + q)jf( P + q y, (23) 

where the color indices are contracted as before. The infrared pile-up of eigenvalues of T>a 
that we have seen in Coulomb gauge must translate into the infrared enhancement of the 
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ghost propagator in Landau gauge. With an implied ensemble average over the gauge field 
the ghost propagator should take the form 

^»-- ( _ p2 _ ze)1+K p 2 «A 2 . (24) 

Here we have reintroduced the scale A. k = 1/2 is a special value since k > 1/2 leads to 
an infrared divergence in IP^g 2 ) as q 2 — > 0. In three space dimensions the enhancement 
corresponding to k = 1/2 gave a confining linear potential and was more than sufficient to 
produce the infrared divergence. But in four space dimensions we would need k > 1/2 to 
produce an infrared divergence in the Coulomb self-energy. 

In section 2 we discussed how the rapidly growing numbers of Gribov copies can cause gluon 
propagation to be highly damped in the infrared, such as G(p 2 ) ~ —p 2 /A 2 for — p 2 <C A 2 . 
In this case gluonic effects can be integrated out on scales of order A, giving rise to effective 
nonrenormalizable interactions at this scale. As discussed in [17] this disrupts a naive power 
counting, where it is assumed that momenta much smaller than A can dominate all loops in 
the infrared. Due to the suppressed gluon propagator, loop momenta of order A can dominate 
some loops. These effects can induce effective couplings of the ghosts to other colored fields 
such as j gh/i j M /A 2 where j ghfl = c a d tl c a and j M = q a, "f fl q a . As well there are self interactions 
among the ghosts, such as the interaction j ghM j gh /A 2 with Ai « A 2 ~ A. (There are also 
interactions among the color octet versions of the currents.) This has the consequence that 
ghost loops, each described by LT"^(g 2 ), may be inserted into diagrams, with each insertion 
contributing an enhancement factor of order (A 2 /(— q 2 )) 2 ^ 1 for k > 1/2 [T7|J§ Due to the 
ghost coupling to the quark current, this can give rise to a power-law enhanced (or log enhanced 
for k = 1/2) anti-screening mechanism. For instance the Wilson loop should be dominated by 
multiple ghost loop contributions. Perhaps this ghost loop anti-screening mechanism could 
be developed further as a covariant description of confinement. 

In conclusion, we have explored the problem of gauge copies present in the continuum 
definition of non-abelian gauge theories and found that an explicit counting of these copies 
leads to insights concerning the infrared behavior of these theories. In the last section we 
found a pile-up of small eigenvalues of the FP operator. Such an enhancement has been widely 
discussed in the literature, and in the Gribov-Zwanziger approach it is said to arise due to 
the constraint on configurations. It is argued that configurations will typically lie close to 
the Gribov boundary where some eigenvalues of small, due simply to the large dimensionality 
of the configuration space [TJ US]. We are instead finding a power-law enhanced density of 
small eigenvalues by studying the numbers and distributions of copies, in the absence of the 
constraint on configurations. 

6 This reference also discusses how ghost loop effects can feed back and affect the gluon propagator. 



16 



Small eigenvalues correspond to small curvatures and thus enhanced gauge fluctuations in 
the direction perpendicular to the gauged fixed hypersurface. In sections 2 and 3 we described 
the infrared suppression of the "physical" gauge modes within the gauged fixed hypersurface 
due an effective shrinking of the configuration space in these directions. This is controlled 
by the rapid growth of N(A). Thus the potentially physical modes are suppressed while 
unphysical modes, modes that cannot appear in physical states, are enhanced. These two 
phenomena together present a picture of confinement. The scale A characterizes the mass 
gap of the theory, the scale below which unphysical modes dominate the dynamics and above 
which the perturbative modes become evident. 

In section 2 we also discussed a vanishing contribution to vacuum energy due to the 
discreteness in the number of Gribov copies. This result is of interest given the present lack 
of evidence for vacuum energy in QCD with massless quarks (THj. We note that the operator 
product expansion in (TI21) also shows that a quark condensate does not mimic the effect of 
a gluon condensate in its effect on the gluon propagator. This is a consistency check on the 
possibility that the nonperturbative physics responsible for chiral symmetry breaking also does 
not contribute to vacuum energy in massless QCD. In this respect it would be the same as 
another possible nonperturbative contribution to vacuum energy, that from instantons, which 
is known to vanish in massless QCD regardless of chiral symmetry breaking (i.e. vacuum 
energy does not depend on the vacuum angle 6 in massless QCD) [20] . 
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